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This report serves to document the LinCom effort under Exhibit B of 
Contract NAS9-lt>097 in supporting the JSC study of the use of the GPS 
navigation system on the Space Shuttle and in Ku-band problem 
investigations. LinCom has been tasked to evaluate system 
implementation, f^rfonnance, and integration aspects of the Shuttle GPS 
system and to provide independent technical assessment of reports 
submitted to JSC regarding integration studies, system studies and 
navigation analyses. The results of this effort is contained primarily 
in a set of five memos [1]-Cb] and a report [bj. In addition, an 
ongoing effort is to review and evaluate the DoO Joint Program Office 
GPS Phase I! User System Specification [7] and the Space Shuttle Orbiter/GPS 
Procuren»nt Specification [8], Documents [1] to 16J are attached for reference. 


^ -0 ef fort REIAT ED TO GPS 

!n s'.ippor* of the Ku-band program, LinCom was tasked to evaluate 
the effect of coaxial cables and connections on the data waveform, to 
perfon.i an analytical investigation of bit synchroniser and phase'Offset 
correction circuitry, and to Investigate discrepancies between predicted 
and measured performance parameters as test data becomes available. 

1.1 Conference Support 

LinCoi’' ['art icipated as a representative of JSC in the GPS 
Const el 1 at i oii Conference at S"diS0 on June S, 19t'>0. The purpose of the 
GPS Constellation Conference was to address the budget-dri ven progr.i-’ 
vh.ange to reduce the number of satellites froin J4 to IG in the Gl'S 
nelworP . Its impact on the Space Shuttle program in terms of Shuttle 
mission schedules and Shuttle navigation via GPS were docuciented in [1]. 


Uc did nut attend the Magnavox Phase II PDR due to attendance space 
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limitations. However, we have reviewed the POR notes (taken by Or. Jim 
Rawlowski and Howard de Vezin) and provided preliminary comments [3]. 
Since the two-channel set Is functionally similar to the Shuttle 0’S 
R/PA baseline, some effort was spent on Interpreting and understanding 
the TTFF procedure for the two-channel set and the results were 
documented. 

1.2 GPS General Support 

Support was provided to JSC on analytical and/or technical 
questions regarding GPS development implementation and performance. 
These efforts are highlighted as follows. 

The timing structure of the received GPS signal was summarized to 
highlight its impact on the normal acquisition procedure and time-to- 
first-fix (TTFF). In particular, the handover procedure from C/A-code 
to P-code was addressed. The results were documented in [2], 

The problems and options associated with testing the GPS receiver 
were investigated in [5], In particular, the emphasis was on testing 
the GPS receiver when only one simulated emitter signal (out of a total 
of a complete set of four) is available. The ideal test configuration 
and the basic system and subsystem tests required to determine receiver 
performance were described. The implications of testing with a nonideal 
configuration, for example, with 1 satellite simsilator, were then 
considered. Special attention was devoted to the 2-channel sequential 
set (present JSC baseline). Supplementary tests that can be performed 
independent of the simulator signals were also highlighted. 

A simple sequential one channel set was recoesnended for JSC 
breadboard developiient. The rationale was documented In [3]. 

In LA], a hardware and software method for obtaining bit and 
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subframe synchronization were proposed. The hardware method was 
reconwended for the GPS breadboard development program. 

1.3 Analytical GPS Modelling 

In order to support the hardware development and Implementation of 
the bit and frame synchronizer [4], we have analytically modeled these 
synchronizers. The results from this analysis a**? provided in [6] and 
should be useful in setting system parameters pertinent to their 
acquisition behavior and detection performance. 

1.4 Record, lendat ions for Additional Areas of Investigative Effort 

LinCom should continue to identify and study key system parameters 

relating to the performance of the Shuttle/GSP R/PA such as: (a) Time- 

to-first-fix (TTFK), (b) acquisition strategies and (c) various tracking 
loop uarameter for PN tracking, carrier tracking, bit synchronization 
and frai.it' syncfironizat ion. One t;iajor purpose of the analytical study is 
to parallel and support the review and critique of the R/PA procureiaent 
specifications as they evolve. LinCon should also continue to follow 

the uPS JPO Phase II efforts in order to keei> track of current develop- 

aients and identify techniques, hardware problems and areas of concern 
which are of interest to the Shuttle GPS program. 

Two additional areas were uncovered during this phase of our 
efforts and they deserve further investigation. One area is the study 
of the i'/PA systcii perforraance testing requirenents and procedures. In 
i articuld'', an u'nk*rstandi mj of the required lest eqiiipments, e.g. CPS 
sK'ioal sii’nil vitO'-s and R/PA test access points are needed. Another area 

is cu'Kernod with the hardware design support for the ISTL GPS receiver 

develop. sent I'rogras' undertaken internally at JSC. 


XinCoi 
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LinCoin has developed an analytical sinwlation of the RU-142 coax 
cable which supports data transmission on the 2-to-50 Mbps channel. The 
analytical sinwlation accounts for cable attenuation, input asymmetry, 
data/clock offset, input amplitude variations, rise and fall tiiros, and 
transition density. The simulation was developed in such a manner that 
the cable pulse output features can be assessed and used in the 
developitent of a specification to be supplied to the tCu-band vendor. A 
time error budget was established in support of the new design of the 
mid-bit data reclocking circuitry. The results of this simulation are 
documented in Section 3 of this report. In addition, LinCom supported 
numerous Ku-band technical reviews at Hughes Aircraft and Rockwell 
International, as required by JSC. 

Z.Z Ku-Band Aiialvtical Investigations 

LinCop! has performed an analytical investigation of the bit 
synchronizer and piiase-offset correction circuitry presently recommended 
by the Ku-hand vendor. Particular attention was given to the false-lock 
problem to which the current design is susceptible. The cable output 
distortions and their impact on the performance of the Ku-band mid- 
bit/phase correction circuitry were evaluated. The results of this 
study are docunented in Section 4 of this report. 

J.3 Ku-Band Test Data evaluation 

Li'ii,nn tias perfumed investigations of di sc^epaficies between 
theoretical pertornance predictions and measi:red values. To this end, 
LinCoii analyzed the available test results, performed analytical 
perfontiance predictions and reported on the major discrepancies and 
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3.0 R6142B CABLE RESPONSE TO SO MBPS DATA STREW IH THE PRESENCE 


Introduction 

An analytical ccmputtr siHKilatlon mcKltl for tfw transfer functlcm 
of a coax cable was develo^d In a previous LInCoro raemo^^^. The 
objective of this report Is to extend the previous model to a case 
where coaxial connectors are present. The analytical model of a 
connector for bit rates under ccmsideratiwi Is developed using 
the manufacturer supplied voltage standing Mve ratio associated with 
tne connector. Tf« connector and the cable ffK)de1$ are employed In a 
software package to study and predict the signal distortions caused 
by the line. It Is sh<wn that presence of a few connectors (<10) will 
contribute very little to the amplitude loss or other distortions 
introduced by the line. 

1. Analytical Model forRS142B Coaxial Transmission Line In Presence 


oaxiai connectors 
For data bit rate of 2>50 ^tt>p$, a coaxial cable can be modeled as 
a linear system with transfer function C(w): 


C(«) ■ e 


where Y is a constant depending on the cable and t is the cable length. 
It Is shown that If (1) Is utilized for GR142C cable, the amnlltude 
loss at 50 MHz equals 2.7 To consider the connector role, one 



can take advantage of the specifications chart prov1<^ by the cwwiectm* 
manufactyrer. At the bit rate under conrl deration, the asotv .leful 
Infomatlon provided for a connector Is the voltage standing wave 
ratio (VSUR) because using this ratio the reflection coefficient can 
be computed. From the tests performed by different nmnufacturers. It 
1$ apparent that the VSUR of coaxial connectors Is linear In frequency, 
and, for the sake of simplicity, one san^le Is selected to be used In 
this report, ^cording to SOCAPEX^^^ catalog a SHA connector for 
R6142B/U cable denx)nst rates the following VSWR: 

VSWR ■ 1.1 ♦ .OZF (2) 

where F Is the frequency In GHz. The relationship between VSUR and 
the reflection coefficient p Is given by 


1p| 


VSWR-1 


(3) 


To derive the transfer function equation associated with a 
connector wit reflection coefficient p, we observe that the reflected 
wave contains a reflected p(wer term given by^^^ 

■IpI^P. M) 


where P^ Is the Input power. Transmitted power P {power passing 
through the connector) can be computed as 

P • % - P. • 0-lol^)P4 (S) 

{41 

Since the power Is proportional to the square of the voltage' , thus 
the folltxvlng equation can be deduced from (S) 
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where V Is the transmitted sMveform voltage and Is the Input voltage 
The ratio V/V^ defines the connector transfer function H(w)t 

H(w) ■ (6) 

Equation {6} Is the flMXtol for a coaxial connector. above i^el 
does not take into account the R.F. leakage or any other losses since 
the frequency of Interest Is rather low (considerably less than 1 GHz). 

Because of the linearity of the overall syston, a cable of length t 
plus n connectors demonstrate the fo1i(»«1ng total transfer function 

• C(«)h"{«) - e‘^^0-|p|^)”^^ (7) 


where |p| Is given by 


|P| 


.1 ♦ .02F(GHz) 


( 8 ) 


Eqjat1(W) (7) Is used In a slimilatlon program to predict the system 
response to different Input patterns. Figure 1 displays the Bwgnltude 
of iij.(w) for two different conditions; plot a Is obtained asswilng no 
connectors on the line, while plot b assumes 9 connectors are present. 
At SO MHz, the amplitude loss Increase because of 9 connectors Is 
approximately .1 dB. 

2- Tot al System Response 

Figure 2 shows the system response to a trapezoidal shape pulse 
of width 20 ns and 3 ns rise time. The three plots correspond to no 
connectors on the line, 9 connectors on the- line, and 99 connectors on 
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FREQOEHCY IN MEGAHERTZ 




Time * 3ns). 
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line. The fir&t tm plots ere very close to tech other irfiich is en 
indication that the distortion intrtkiuced hy 9 connectors is negligihle.^ 
Figure 3 displays the output to a square Mavefora with 3 ns trapezoidal 
rise tin». Here, only the two cases of 0 and 9 connectors are plotted 
and the two plots are almost indistinguishable. Figure 4 illustrates 
the system response to a random pattern. Again the plot for 9 connectors 
is indistinguishable from the one with 0 connectors. The plot for 99 
connectors indicates a significant Increase in signal distortion. Figure 
5 illustrates the output rise time versus the input rise time. Plot a 
is obtained using the 30-70% definition while plot b corresponds to the 
20-80% definition. Using plot b, 5 ns output rise time is measured when 
the input rise time is 3.5 ns, while if there were no connectors, the 
input rise time could be as high as 3.8 ns^^^. This means that the 
appearance of 9 connectors on the line restricts the input rise time by 
an extra amount of .3 ns which, of course, is small. However, in case 
of 10-90% definition the change in the input rise time is greater than 
1 ns. This last observation is one of the many reasons why 10-90% spec 
is not suitable for rise time definition. 

3. Amplitude Loss 

Figure 6 displays the amplitude loss versus the signal asymmetr-. 

Nine connectors are assumed to be present. If Figure 6 is compared with a 
similar figure in Reference (1), about 1% increase in the amplitude loss 
becomes apparent. Plot a of Figure 7 displays the line response to a 
20 ns pulse while 9 connectors are present. Plot b shows the output 
if 25% asymmetry is imposed on the pulse. The input rise time is 3 ns. 
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Figure 3. Square Wave (Input Rise Time 3 ns). 
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Figure 7. 
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4. Threshold Variation 

Figure 8 displays the timing offset generated by the threshold 
variations. After cwnparing this figure to a similar one In Inference 
(1)» the fo11(Mlng observation Is made. The timing offset Is not 
Increased by the presence of 9 connectors whm threshold variation Is 
less than 10%. For a variation between 10 and 25%, the offset Is 
Increased by a small amount; this Increase, however. Is not large 
enough to be of any practical Importance (less than .03 ns). 

5. Conclusion 

This report Is a continuation of the work done In Reference (1). 

The model utilized to represent 100 ft of R6142B cable introduces 2.7 dB 
amplitude loss at 50 MHz. It Is shown that the addition of 9 connectors 
to the cable at 50 Mbps data rate Increases the output degradation by a 
negligible amount. The threshold variations cause an unwanted time offset 
at the data reconstruction mechanism output. The following table shows 
the timing error for two extreme cases. The first case corresponds to 
a source putting out pulses with .5 volts as the lower level and 6.5 volts 
as the higher level. The second case coincides with a source putting 
out -.5 volts and 4.5 volts as two pulse levels. The normal situation is 
considered to be 0 and 5 volts for the two pulse levels with the threshold 
set at 2.5 volts. The table also contains the amplitude losses related 
to the sources under consideration. 

Note that the threshold variation Is computed In percent so that 
Figure 8 can always be used. For example, to compute the threshold variation 
for source 1, number 2.5 Is subtracted from 3 (source average) and Is 
divided by the high (6.5) and low (.5) voltage difference (6), and the 
result is multiplied by 1(K). Amplitude loss is also represented in 
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percent In Figure 6. To convert the amplitude loss to volts the nunfter 
(Stained from Figure 6 Is divided by 100 ai^ multiplied by the high and 
low voltage difference. 


Table 1. Source Amplitude Variation (3 ns Input Rise Time). 


SOURCE 

SOURCE 
LEVELS 
IN VOLTS 

THRESHOLD 

VARIATION 

TIMING 

OFFSET 

OUTPUT AMPLITUDE LOSS 
(20% ASYMMETRY) 

1 

.5-6.5 

Vt 

♦.9 ns 

.75 volts 

2 

(-.5)-4.5 

10X 

-.5 ns 

.625 volts 
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Introduction 

The uncocted 2 to 50 Mbps data must be corrected for asymmetry 
before being processed by the convolutional encoder In the Ku>band 
return link communication subsystem of Space Shuttle Ortiter [1]. 

This report describes the function of the asymmetry corrector hard«vare 
designed by Hughes and also verifies Its timing budget. The corrector 
circuitry Is under construction at the present time and Is approaching 
Its final phase. 

1. Timing Allocation 

A total of 5.9 ns Is allocated to three iimjor causes of asymmetry: 

a) Source asymmetry T, defined as the difference between a 1 pulse 
and a 0 pulse width. 

b) Data rise and fall time contrllHitlon due to amplitude loss T|^ 
defined as (20-80% definition for rise time): 

T « (^plltude Loss )x (Rise Time) 

'b lOOx.6 

c) Sampling margin reduction due to error In amplitude threshold 
estimation T^. Figure 1, copied from Reference [2], can be used 
to ccKnpute this tin«; for 13% amplitude loss T^ equals 1 ns. 

The three timing factors are computed for the worst case of 50 Mbps 
data 20 ns pulse width. The rise and fall time definition Is 20-80% and 
the rise (and fall) time spec Is 6 ns. Amplitude loss (because of the 
cable) 1$ 13%. 


Using the above definitions, time T. Is computed first. 
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Figurt 1 can bt used to determine r 

€ 

. ins 

To determine t)M maximim alloMable source asymnetry A, mist be 
com^ted first: 

Ta - 5.9 - T|, - Tg ■ 3.6 ns 
Since tte (Hilse width Is 20 ns» 

A ■ ^ X 100 * 18S 


For safety reasons A Is selected to to 16X. 

Table 1 sunnarlzes all the contributing factors to loss In san^llng 
margin, /^cording to this table the worst case (roaxlnnmi) loss Is 9 ns 
or 45t ("2j y— ). The mid-bit corrector harAnare functions properly only 
If loss Is be1<n« 50%. 

To correct for asymmetry a pair of complementary and synnwtrlc clocks, 
I and Q, are generated, and the data Is san^led by the I clMk at mid-bit 
points while the Q clock Is made to coincide with the rising edges of the 
data. Figure 2 Illustrates the sampling methto on a *10110** data sequence 
with 40% asymn^try (loss In sampling margin). Tto ^thod will work only 
If this asymmetry Is bel<nr 50%. The concept Illustrated In Figure 2 Is 
utilized In asymmetry corrector circuitry. The worst case sampling 
margin Is 5% or 1 ns. 

2. HId-BIt Corrector 

The purpose of the Mid-Bit Corrector circuitry of Figure 3 Is to 
generate two ccxnplementary clocks, called I and Q, and to position th^ 
so that the leading edges of the 0 clock coincide with fN2s1t1ve-go1ng 
data transitions, and the leading edges of the I clock can be used to 
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reliably sample Incoming data. This process Is accomplished by Injecting 
a slowly varying dc volta^ Into the existing phase-lock-loop's voltage 
control amplifier to adjust the phase of the VCO output relative to data 
transitions. 

The phase-correcting volta^ Is generated by a six-bit dlgltal-to- 
analog (D/A) converter. This allows the clock to be phased with respect 
to the data with a precision of one part in 64 (six-bit resolution). This 
precision Is equivalent to 5.6 degrees (360/64). 

The principal components of the corrector circuitry are the phase- 
lock-loop components (phase-frequency detector, amplifier, voltage- 
controlled oscillator (VCO), and dIvIde-by-two flip-flop); the six-bit 
D/A converter; a modulo-64 up/down counter; and logic to detect clock/data 
skew. 


The voltage level out of the D/A converter is controlled by the 
state of the modulo-64 up/down counter. The counter is clocked by a 
1.0 KHz clock, but it is prevented from being incremented at more than 
a 2.0 millisecond rate by external logic. This timing, shown in Figure 4, 
allows the phase-lock-loop to settle from a D/A voltage change for one 
millisecond before the next clock/data phase comparison is permitted. 

The logic works as follows: 

The counter will be Incremented by alternate 1.0 KHz clocks if a 
positive-going data transition has occurred during the previous one 
millisecond time period, which Is indicated by the true state of the 
counter enable flip-flop. During the one millisecond period after each 
counter update, the counter enable flip-flop is held reset by an 
exclusive- OR gate that detects that the counter's leads significant 
bit (LSB) has just changed. As shown, the LSB output is delayed for one 
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Figure 4. Modu1o-G4 Up/Domrn Counter Timing. 
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millisecond through the LSB delayed FF, the output of which is exclusively 
ORed with the LSB output to produce the loop-setting time interval. 

Figure 5 shows the state diagram of the counter enable flip-flop. States 
0 and E are unstable. Starting from state A» upon arrival of a data 
transition, a switch to state 0 occurs. Since D is unstable, the flip-flop 
autcxnatlcally moves to state C. State C is terminated v^en a 1 KHz clock 
transmission occurs. The next state B lasts 1 ms until another 1 KHz 
clock transmission (positive) arrives. This will put the flip-flop back 
into the original state A. While at state B, it is conceivable that 
because of occurrance of data transitions this state would temporarily 
switch onto state E. However, since E is unstable it will immediately 
change to state B. This phenomenon can also happen between states C 
and 0. 

The flip-flop which divides the VCO output frequency by two produces 
complementary I and Q clock outputs, as shown in Figure 6. Incoming data 
is clocked into flip-flop FDI by the mid-bit I clock, then delayed a half- 
bit time through FEQ, clocked by the Q clock. Whenever a data transition 
occurs, the two flip-flops will become out of phase for half a bit time. 

An AND gate detects positive-going transitions which are clocked into the 
counter enable flip-flop by the Q clock provided the reset pulse generated 
during the loop-setting inte*-vai is not present. Thus, the first positive 
transition after the reset pulse disappears will cause the counter enable 
flip-flop to be set. 

To determine if the phase of the VCO clock should be advanced or 
retarded, data is clocked into flip-flop FDQ by the 0 clock. If transi- 
tions occur early (solid lines in Figure 6), FDQ will change a half bit 
time before FDI. If transitions are late (dotted lines in Figure 6), 
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Figure 6. Mid-Bit Corrector Timing. 
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FOQ will change a half bit time after FDI. This phasing of FDQ and 
FOI Is compared In the Up/down AND gate. A logic 2ero at Q clock time 
calls for a clock phase advance; a logic one* a phase delay. 

The sU.e of the up/down gate at Q clock time Is monitored by 
flip-flop Ff.i . When a positive transition Is detected by the counter 
enable flip-flop, the latest up/down gate Information will be stored In 
FHQ and preserved against subsequent changes by the counter enable 
flip-flop outrut state v/hich changes at Q clock t1ii«. The true state of 
the counter enable flip-flop allows the up/down counter to be Incremented 
by the next 1.0 KMz clock; the state of FHQ controls whether the count 
Is advanced or retarded. 

After the loop achieves the correct clock/data phasing, the counter 
will alternate between an up convnand and a down command on successive 
updates. When a new set of clock and data channels Is selected, a 
maximum of 31 new counts could conceivably be required by the counter 
before proper phasing Is achieved. If data transitions occur as often 
as two In any one-millisecond time period, the loop will acquire proper 
phasing within 62 millisecond when 31 counts a.*e required. 

Data sampled by flip-flop FDI will be used by the Convolutional 
Encoder logic. For mid-bit sampling of this data, the Q clock will be 
supplied, along with the X2 clock. 

To reduce power and parts count, only that portion of the mid-bit 
corrector circuitry that must operate at megahertz frequencies has been 
implemented by ECL parts. The low-speed up/down counter, for example, 
will be mechanized by low-power TTL circuits, with level shifters being 
used between different logic levels. 
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3. Threshold Estimation and the Data Receiver 


The purpose of the circuit at the upper left corner of Figure 7 1$ to 
reconstruct sharp bit pulses from degraded Input waveform. Figure 8 
displays a PN sequence after traveling through 100 ft of coaxial cable. 

To recover the original pulse shape, It Is necessary to establish a 
threshold to differentiate between the zeros and the ones of the 
sequence. To generate the proper threshold, the circuit detects the 
high and the low peaks of the waveform and forms the average of these 
two quantities. The .1 F capacitor preceded by a diode responds 
quickly to the proper peak and obtains the peak value. The buffer 
saves the peak voltage for at least 1 ms. When a new source Is 
selected the high peak might drop as much as 2 volts (Table 2). The 
circuit Is relatively slow in responding to sudden high peak drops; 
according to Mr. Pfiffner of Hughes It takes the circuit 350 ms to 
respond to a drop of 2 volts. This observation Is Illustrated In 
Figure 9. 

Assuming an RC model, the time constant of this circuit can be 
computed as follows. Again, observing that 350 ms Is needed to drop from 
6.5 to 4.5 volts and calling the time constant t 


6.5 

O 


e 


.35/t 


or T» 1. With such a high time constant, the voltage change In 1 ms 
Is negligible. 

4. False Lock 

Jn the transient period of frequency acquisition, the phase-locked 
loop tracking the Input clock In the asymnetry corrector may false lock 
if the dc voltage injected into the voltage-controlled amplifier Is 
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not tquti to zero. To ovoid foist lock* the circuit ot tht lower left 
comer of Figure 7 Is used to reset the up/down counter to zero every time 
the clock frequency chonges. The circuit operotes on this principle thot 
ot the end of every 2S6 counts If iv of the six most significant bits of 
the modulo 256 counters Is not o one the up/down counter Is outomoticolly 
reset. This condition con only hoppen If the two frequencies of the 
Input clock ond o clock ore at least 1 ,67%(100/64) apart. For example. 

If Q clock frequency Is 40 MHz ond the Input frequency Is larger than 
40.7W1Z this circuitry will reset the up/down counter. For frequency 
chants smaller than 1.6X, the |^ase*1ocked loop automatically corrects 
the i^ase, hence there Is no need to reset the up/down counter. The 
time required for the circuit to detect a not-ln-lock situation can be 
computed. It is clear that 4096 pulses are required before the up/demm 
counter can be reset. In the vicinity of 2 MHz, this time equals 2 ns; 
while at 50 MHz it is much shorter and equals .08 ms. 

Re q ui sition 
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The acquiisiton process is basically composed of two parts: frequency 
and phase acquisitions. Figure 10a Illustrates the block diagram of the 
phased-locked loop. The phase-frequency detector can operate either in 
a frequency detection (when frequencie;:. are unequal) or in a phase 
detection (when frequencies are equal) mode. Figure 10b shows the approximate 
frequency detector characteristic while Figure l(k illustrates the phase 
detector characteristic. Frequency detector characteristic is only valid 
if frequencies f^ and f^/N of Figure 10a are not equal. 

5-^ Frequency Acquisition Time 

The differential equation describing the loop behavior in the 
frequency acquisitions mode is 
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Figure 10a. Phase-Locked Loop. 
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Figure 10b. Frequency Detector Characteristic. 
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Figure 10c. Phase Detector Characteristic. 


PHASE DIFFERENCE 
(RADIANS) 






where is the (nonzero) phase error - e^/N. Kp is « constant in 
volts that depends on the frequency detector, Ky is the VCO gain, T 
is a constant parameter of the loop filter, and N is equal to 2 . Assum 
ing the initial phase is 3 q and initial frequency offset isAu/2w, the 
first (or second) equaion in (1) can be used to determine (t). 


cf(t) 


= 6fj + (&U, - “ 3T7 ° 


2 


(2 


The frequency acquisition time can be computed by setting the 


derivative of <f(t) equal to zero. 


t / - Ao) • — T . 


TT 


2T, acq 


Thus 


KpKyTi 

■ irr~ 


acq 






Since KpKyT^/T^ is small, equation (4) can be simplified to 

T .Aw 


acq KpKy 

rcr 


(3 


(4 


(5 


Consulting Motorola MC1204C and MC1658 data sheets, typical values 


of Kr and K., can be found. 
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.35 volts 

45x2 xlO® rad/s/v 
16.5x10^° 

“2 

Figure 11 illustrates*^ (t) for maximum frequency offset of 50 MHz. 

^acq computed either using equation (5) or Figure 11. The 

ccxnputed value for is approximately 1.9 ms. Thus, frequency 
acquisition occurs in a relatively short time. 

5.2 Phase Acquisition Time 

The linear characteristic function of the phase detector is illustrated 
in Figure 10c, and the linear loop equation is 




Wl 


+ 




(7) 


where K^is the phase detector gain in volts/rad. Defining w 

T 

and t “ equation (?) simplifies to 

cp + 2w^c = 0 (8) 

Assuming « .111 v/rad and using the data given earlier in this section, 
01 ^ and t can be ccxnputed: 

* 2.29x10® 

C * 3.44 




Solving the second order differential equation(6): 



Frequency (W2tt); £o) = ZwxSOxlO 
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. r ^-.344txl0® . ^ -15.4txlO® 
«f(t) * c^e ♦ c«e 


where end C 2 depend on the Initial conditions: 


Cf(0) « 6/ 


^0) * Au - 2 u|^c6o 


Figure 12 Illustrates <|(t) as a function of u^t. From this figure one can 


conclude that phase Is acquired in less than 1 ms. For example If Au - 0, 


the phase <|f(t) beomces less than .5 degree In just .1 ms. 


After the clock frequency and phase are acquired, the modulo 64 


counter will function as a slow phase corrector to synchronize the Q 


clock with the data. The counter assumes one out of 64 possible states. 
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Frequency Step of 
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since the threshold estimation end other acquisition operations can 
happen simultaneously. It seems a total of 350 ms would be adequate. 
Yet, for safety reasons, a total acquisition period of 400 ms Is 
suggested. 
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GPS-80-CMC-2 

LINCOH MEMO 

To: Don Pusch Date: June 9, 1980 

From: C. M. Chie and W. C. Lindsey 

Subject: GPS Constellation Conference at SAMSO, June 5, 1980 
cc: J. MacLeod, B. Batson, J. Johnson, W. Teasdale 


Abstract 

The purpose of the GPS Constellation Conference is to address the 
budget-driven program change to reduce the number of satellites from 
24 to 18 in the GPS network. Its impact on the Space Shuttle program 
in terms of Shuttle mission schedules and Shuttle navigation via GPS 
are summarized. 

Current GPS Status 

The conversion of the GPS program from concept validation to full 
scale engineering development was approved by DESARK II in 24 June 1979. 

Two user set contracts were awarded to Magna vox and Rockwell Collins from 
which one will be selected. Another contract was awarded to Rockwell 
International last summer for satellite replenishment. The control segment 
contract will be awarded in September of this year. Currently there 
are six space vehicles (SV) on orbit and 4 of them are functionino up to 
specification while the two earlier ones have clock problems. 

Last Christmas, a reduction from 24 to 18 5V in the GPS network 
was directed by OSD as a result of budget cuts. The decision from 
24 to 18 is strictly a cost-saving measure and OSD was aware that a 
constellation of 18 SV is neither optimal nor desirable to support 
civilian and military requirements. 
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Shuttle Launch Schedules 

All 18 SV will be launched by the Space Shuttle either durinn 
shared and/or dedicated missions. This will occur durinq the time-frame 
from late 1984/early 1985 to 1987. In 1987, GPS will be fully operational. 
Each launch will carry 2 to 3 satellites depending on the constellation 
selected and whether the Shuttle Is shared. 

Shuttle Navigation via GPS 

The plan now Is to keep 5 functional SV until 1984 when the first 
Shuttle launch Is scheduled. The new SV will be added In such a 
fashion so as to optimize world-wide coverage percentage. The ultimate 
goal is to select a satellite constellation that allows and maximizes 
3-0 continuous world-wide coverage when GPS is fully operational. A 
degradation of performance in the form of isolated "outage" pockets 
where the navigation accuracy is unacceptable is anticipated from the 
reduction of the nutiHjer of satellites (see next section for more detail). 
For a number of users the exact location of these pockets is more 
important than the total percentage of world-wide coverage. For example, 
the Space Shuttle may need GPS navigation most durinq deployment and re- 
entry near Florida, California and Texas. As a result, the Shuttle program 
may want to have an input to the determination process of placino 
these pockets or to change the criterion of optimization in selectina 
the satellite constellation. 

Impact on Navigati on Accu racy 

The 18 SV constellation is expected to provide degraded ner<o»'"!ance 
relative to the original 24 SV network proposed. Durinq certain ti.,e 
periods in isolated geographical regions, outages of GPS receiver's are 
experienced. By outage is meant that the 3-D 1 range error exceeos 7D 
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meters assuming a 1-D l<j range error of 7 meters (P-code used) and 
10 degree elevation mask angle. The nominal In value Is 16 meters. In 
the 3 orbit x b SV constellation studied by Aerospace Corporation, these 
outages last for about one hour 3 times a day In four major Isolated 
areas on the oceans. Each area Is a little over the size of CONUS. 

If an altimeter Is available with 1o accuracy of 30 mters, the outage 
areas can be drastically reduced due to an improvement In nwasurement 
accuracy. 

As an example, the following chart was shown: 


Altitude Error 

Accuracy 


3D 

2D 

300M 

335M 

148M 

30M 

36M 

21M 


The altitude error Is the 1o error of an altimeter. With an error of 
30W‘, the 3-d1n«ns1ond1 position error Is 335 meter and the 2-dimensional 
posit'ion is 148 meter for the worst case outage region. The accuracy Is 
reduced significantly If a better altimeter Is used. 

Two other constellations are also studied: 18 orbit “rosette’ 

(Phllco Ford) and 6 orb1tx3 SV (STJ). Similar outage patterns can be 
observed although the areas affected are much smaller and scattered. It 
appears that the 6 orbit x 3 SV constellation offers the best performc .ce 
found so far. As a simple comptrlson the following chart was shown: 
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Percentage of Time Total Number of Satellites Visible 


Constellation 
(Number of 
Orbital Planes) 

5 Itegree Elevation 
Angle Mask 

Number of Satellites 

10 Degree Elevation 
Angle Mask 

Number of Satellites 


4 

5 

6 

7 

4 

5 

6 

7 

3 

100 

92 

72 

40 

98 

7 

7 

7 

18 

100 

77.6 

77.8 

30.8 

99.9 

86 

44 

7.6 

6 

100 

99 

78 

28.7 

100 

90.4 

41 

5.9 


Since 4 satellites are required for a normal position fix, the 6x3 
configuration appears to achieve the best worst-case performance with 
a 10 degree elevation mask. 

At this point, the ongoing studies for optimal constellation 
selection are driven by minimizing the percentage of outage areas 
on a world-wide scale. It appears that by nwdifying the orbits, these 
outage areas can be moved around. 

Unctassified charts presented will be available to LifiCom shortly 
and any pertinent data will be forward to JSC as soon as possible. 
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LinCom Corporation 


Date: Noveif^r 4, 1980 


From: C. M. Chle 


cc: Bart Batson, Don Pusch, 
Jack Johnson, Ann Sullivan 


To; Jim PawlOMSkI 

Subjact: OPS Timing Structure and Typical Acquisition Sequence 


Stmitiary 

This memo simwarlzes the timinq structure of the received GPS 
signal and Is simplified to highlight Its impact on the normal 
acquisition procedure and t1me-to-f1rst*f1x (TTFF). In particular, 
the handover procedure from C/A-code to P-code Is addressed. 

1 . 0 GPS Signal Structure for Time Peference 

A simplified GPS signal timing diagram is shown in Figure 1. All 
the information about the particular broadcasting satellite and the 
almarac of the GPS constellation are packaged in a 30 sec frame of 
data. The data frame is repeated every 30 secs starting from the 
end/start of the i««ek (midnight Saturday night - Sunday morninq). The 
data fran« is divided into 5 subfranies of 6 seconds each. In the beginning 
of each subfrrme there are a telemetry word (TLM) and a handover word 
(HOW) to help facilitate the transfer from C/A code to P-code. The TIM 
data slot contains a preanible and a telemetry message The prea^ible is 


I 

I 

I 

I 


used to resolve the phase ambiguity involved with supnressed carrier 
tracking of the navigation message data which is NRZ coded. The HOW 
slot contains the truncated Z count which is the ti'ie {ti’’ie of wee^ ' of 
the leading edge of the next subframe. The sub'rame ID is used to 
indicate tne kind of information contained in the data block. 


The 50 bps daU Is modulo 2 added (XOR) to the P*code (I carrier 
channel) and to the C/A code (Q carrier channel). The P>code Is a 
composite code cmsisting of and Xg. The period of x^ Is 15,345,000 
chips long and spans 1.5 second. The all "l"'s epoch is llnea un 
Mith the leading edge of the data subframe as well as the C/A code. 

The C/A code Is 1023 chips Iwtg and spans 1 msec. 

2.0 Normal Acquisition Sequence 

In the normal acquisition sequence, the receiver first tracks the 
C/A code. Using the code epoch, It can determine the emitter's time 
to within an ambiguity of 1 msec. The receiver then proceeds to 
demodulate the data. Two pieces of Information are sought at this 
point: the Z -count and tne epoch of the next subfranw leading edge. 

This Information Is used to set the local P-code generator so that It 
can be started at the same code phase as the Incoming P-code at the 
upcoming subframe leading edge. This Is possible since the epoch of 
the P-code, the C/A code, and the leading edge of the subframe all 
occur at the same Instant, Note that if the receiver's clock Is accurate 
to 1.5 secs, the Z -count can be generated Internally by the receiver. 

A detailed acquisition sequence is typically as follows: 

1 . C/ A Code loop Acquis ition and Track 

a) Start C/A code search. 

b) C/A code loop pull in and track. 

J J • D ata Dgnodulation 

c) AFC pull In. 

d) Costas loop pul 1 in. 

e) Bit sync pull in and track. 

f) Ambiguity resolution and subframe sync (waxlmuir l 
subfranie = 6 sec after previous step). 



.. 4JUeiiU.m I1.J..I. .... 






III. P-Code Handover 

g) Obtain Z count and align local P-code phase. Continue 

to demodulate data (need subframe 1, 2» and 3 to obtain clock 
error model, satellite ephemeris; need subframe 4 and 5 
for message and alamanac). 

h) P-code acquisition and track. 

1) Pseudo range and delta range measurement and Initiates 
navigation software. 

3.0 Notes 

The range delay between a GPS emitter to a ground user and a 
second GPS emitter to the san« user can be as high as 0.02 sec. The 
arrival epochs of the signals from three different emitters is depicted 
in Figure 2. Hence the HOW from one emitter cannot be used to pinpoint 
the 1.5 sec P-code epoch of another ©nitter which is lined up with the 
leading edge of the subframe of that particular emitter. 
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Figure 2. Arriving Epochs from Different Emitters. 
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PROPOSED SPEC CHANGE 



3. 2. 1.2. 3.1 Nornial Acquisit io n/State 1 . In this state the R/PA 
shall acquire and track the NAVSTAR C/A signal aided by range and 
doppler estimates computed in the R/PA. When the handover word (HOW) 

IS recovered from the 50 bps data stream the R/PA shall acquire and 
track the P signal. The R/PA shall track the carrier and demodulate 
data from the carrier. Pseudo and delta range measurements shall be 
made to full precision. 

3. 2. 1.2. 3. 2 Direct Acquisition/State 2 . In this state, the R/PA 
shall acquire the P signal without first acquiring its associated C/A 
signal. P code phase and frequency estimates for the acquisition shall 
be derived from almanac or current ephemeris data, present position 
and velocity estimates, and a precise input of Coordinated Universal 
Time (UTC). Direct P code acquisition shall be used to reacquire 
transmitters due to constellation revision and to improve anti-jaimiing 
margin when necessary. 
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INTEROFFICE MEMORANDUM 
GrS-81-CMC-l 

Date: January 19. 1981 
From: C. M. Chle 


To: Dr. Jim Pawlov^ski 

Subi**ct: Magnavox PDR Review and NTC Highlights 

cc: John MacLeod, Jack Johnson, William Teasdale, Ann Sullivan 


SUMMARY 

This memo serves to document a preliminary review of the Magnavox 
GPS Phase II PDR. Since the two-channel set is functionally similar to 
the Shuttle GPS R/PA baseline, some effort was spent on interpreting and 
understanding the TTFF procedure for the two-channel set and the results 
are documented. In a forthcoming memo, we plan to compare the 
difference in TTFF definitions used by Phase II set and the Shuttle GPS 
R/PA Spec. 

During NTC, a paper was presented by Hewlett Packard (HP) for a low 
cost receiver. Judging froii its schedule (approximately one year) and 
simple mechanization, it appears to be a viable candidate for the JSC 
GPS breadboard development program. Even though the receiver is rather 
unsophisticated, it does provides an opportunity to familiarize oneself 
with the operational aspects of using the GPS navigation signals. 
Preliminary Review of the Magnavox GPS Phase II PDR 

Based on the PDR notes (taken by Dr. Jim Pawlowski and Howard de 
Vezin), here are some coments of an overview nature: 

• It appears that the design is geared towards operation in a 
high jaricTiing environment. 

t The concept of modularity/conmonal ity is being stressed in 
order to tailor the set configuration to various host vehicle 
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requirements by simply assembling the receiver from common modules. 

This concept also permits major LSI development for common hardware 
components such as User Time Clock (UTC) and coders, etc. 

• The significant improve<nents on the Phase II sets are on the 
operational aspects rather than functional aspects. For example, the 
following functions are highlighted: 

• Automatic Fault Isolation (AFl) 

• Alert and Status Display 
t Self Test Capabilities 

• In terms of receive*- hardware realization, the five channel 
set is basically an X set; the one- and two-channel sets are based on 
the 2-set/Manpack concept. 

• More attention is devoted to the interface requriements 
between the receiver and the host vehicle system for navigation aiding. 

t It appears that the five-channel set can be adapted easily to 
meet the Shuttle GPS R/PA specifications in terms of dynamics and TTFF. 
Dual Channel .Acquisition Sequence 

The nomal acquisition sequence for a dual channel set in motion is 
depicted in Figure 1. The first half of the acquisiton sequence is 
devoted to acquiring the C/A code. The second half operates on the P- 
code. After initialization, the first satellite signal is acquired 
using coordinated code search with the two available channels. Since 
the tii.ie uncertainty is specified to be 20S (lo), the sequential code 
search is performed uniformly over the 1023 C/A code chips. Once the 
code signal is acquired in one channel (say channel 1 as in the figure), 
the receiver continues to perform: 

• Code Loop Pull-In 




t AFC to Center the Costas Loop 
t Costas Loop Pull-In 

• Bit Synchronization to Line Up the Bit Edge 
t Resolving the Ambiguity Inherent in the Bit Sync from 
the TLM Word 

The receiver then proceeds to demodulate the HOW and sets the user time 
using the HOW. (Notice that the GPS time can be determined from the HOW 
to within the range uncertainty between the user and the satellite which 
is in the order of 35 Km (lo) or 120 C/A chips. Once the user time is 
set, the receiver can then change to a Gaussian 3o search to acquire the 
remaining satellites. The Gaussian search is much more efficient time- 
wise and amount to substantial TTFF savings.) Channel 1 then continues 
to ♦'ake pseudo range and delta range measurements. It then filters 
these measurements and prepositions the VCO for the code signal from the 
third satellite. 

Once the first channel verified the presence of the code signal, 
the receiver on the other channel {Channel 2) starts the uniform search 
for the second satellite. If the first channel can read the HOW and set 
the user time before the acquisition of the second satellite in Channel 
2, the second channel will switch to the Gaussian search node. The 
receivers continue to acquire and measure the rest of the sigr.als as 
indicated in the figure. Before switching tu tracking the P-code. the 
set reacquires the satellite signals again to compensate for any 
untracked motions between the user and the satellites during the gap 
between the channel sequence. After the C/A to P code handover is 
completed, the set proceeds to calibrate out the differential hardware 
delays between the two channels (including Li and frequencies). This 





can be done by tracking the same satellite signal using both channels. 
The set then proceeds to demodulate ephemerides data from the satellite 
signals. First of all, it must wait for the frame sync and starts 
collecting data from the first two satellites during subframes 1, 2, 3 
(18 sec). During su‘ 'rames 4, 5, the set reacquires the remaining 
satellite signals ^nd collects data when the next frame sync comes up. 
The set then makes the ionospheric measurement and starts the sequential 
tracking (or dwell). After one cycle of the sequential track, the first 
fix is available. 

The total Time-to-First-Fix TTFF is 203.1 sec for a set in 
motion. Notice that 60 sec is devoted to ephemerides collectio'i, a 
maximum of 30 sec can be wasted for wailing for frame sync (FS), and 

there is a maximsjn of 6 sec to wait for the TLM word to resolve the bit 

sync ambiguity for each signal. 

Low Cost GPS Receiver by H !* 

The low cost GPS receiver built and tested by HP was designed to 
meet the requirements of land navigation with limited mobility (ie. , 
s^^ationary or up lo GO n^ih). It was initiated in 1977 and built and 
tested successfully in 1978. Apparently, two engineers were Involved 
(level of support?). Tho fiP receiver is a sequential one-channel set 
using C/A code on the Li carrier. The tost results showed sn accuracy 

of 20 meters. A Z-80 nicroco-iputer is used ^or receiver control wb'Me 

the navigatioii computations are per^or icd on an !iP982G dest tor 
calculator (see Fig. 2), The navigation solution for the position .e’d 
cIocL bias is computed itei^ati vely. No atter.pt is made t>> smoo.h out 
consecutively fixes with Kalman fi'terlng. Hence, the computat ions 
involved are basically for solving a four-dimensional nonlinear 


equation. 

Reference 


K. P. Y1u, et al., “Land Navigation with a Low Cost GPS Receiver," NTC , 
December 1980, Houston, Texas, pp. 55,3.1-55.3.7. 
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R1 - Reacquire Satellite* #1 
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INTEROFFICE MEMORANDUM 
RPS-81-CMC-2 

February 17. 1981 
From: C. M. Chle 


To: Dr. Jim PawlOMSkI 


Subject' Bit Synchronization and Message Demodulation for GPS Signal 

cc: John MacLeod, Jack Johnson, U1 111am Teasdale, Jack Seyl. Don Pusch. D. 


S urinary 

A hardware and a software method for obtaining hit and subframe 
synchrcniza*. Ion are described In this neno. It appears that the 
hardware method -’s more appropriate for the GPS breadboard development 
program. 

GPS 'lessane Tinlnn Structure 


Figure 1 shows the timing structure of the received uPS ncssage. 
The message is structured to facilitate data derH.w laMon. The 
Inpo-tant feature Is that the C/A code epochs (and P-rode epochs) are 
coherent with the data bit boundary and each data bit Is equal to 20 
complete C/A code periods. Each subframe starts with a preamble word 
(IJOOlOll) whose purpose Is two-fold: (a) to signal the start of a 
subframe and (b) to provide a means to check the polarity of the 
riemod'i’ a' ed data. 

In this 'icM'j, we assurie the code loop Is In lock, I.e. . the C/A 
'o!': o,v)ch« ir.« »7i1 Table. Before vve proceed to de;io(1ulate trie 
'n/ioa' )i 'fC’ly, we rtist resolve the following a'^’M^ulty 
■ on ii t ’ O’lt 

































(a) code epoch anbigulty - determine which one of the 20 code 
epochs Is coherent with the bit clock. 

(b) data bit Inversion ambiguity - determine whether the 
demodulated bit has been inverted. 

(c) Subframe ambiguity - determine which bit Is the start of a 
Subframe. 

The first ambiguity condition Is resolved by means of a bit sync circuit 
(or procedure). The second and third ambiguity conditions are resolved 
by the preamble word. 

Typical Bit and Subframe Synchronization - Hardware Approach 

A hardware approach for bit and subframe synchronization procedure 
for resolving the ai.ibigulty conditions described In the last paragraph 
is exemplified by the techniques used in the Stanford Telecommunications 
Inc. (STI) Geodetic Systen Receiver. Figure 2 shows the circuit for C/A 
hi* sync and d.ta detection and Figure 3 shows the subframe sync 
procedure. 

bit sync circuit operates by checking for data transitions 
arou.id each epoch (2C possible epochs per bit). There ere 20 
accumulators keeping track of the transitions scored by a particular 
epoch. After one of the accumulators fills up to 2 ^ * 64 counts, the 
asscciatcd ci’och is declared to be the data hit boundary and the data 
cloc* is set arcordingly. The data dock is used to provide the timing 
frr the 1”) 1 1 1 second integrate and djiip circuit for the data detector. 

!n order to simplify hardware, the possible data transiton is 
le'eruMcd ljo*ing at portions of 1/20 bit as shov/n. The transition 
detector shof/n in Figure 2 can be a simple comparator, f.'otice that 
since only 1/20 of a bit is integrated, the signal-lo-noise ratio 
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degrades by 13 dF.. To conpensate for this degradation* the accumulators 
are used. The slgnal-to-nolse ratio Improvement Is proportional to the 
minber of transitions observed before declaring "lock". Further 
analysis Is necessary to determine the appropriate accumulator size for 
a particular operational requirement. 

In Figure 3, the navigation data from the data detector output of 
the previous figure is resolved by matching the data stream with the 
prejrihle pattern. Once a match is found the data 1s Inverted 
accontingly and suhfraine sync Is declared. The HOW Is then demodulated 
and used to set the local receiver time. 

Typical Pit and Subframe Synchronization » Software Approach 

A software approach for bit and subframe sync is typified by the 
Magnavox X-set. the principle of operation Is very similar to the 
hjrdv/are apprudch. Cxcerpts enclosed describing the procedure are from 
an Air F jrce Perort SAfISO TR77-102 Vol 1 entitled A Functional 
jcscription of the Navstar GPS Recower ttodcl X. No»e that the software 
aiiproach is appropriate here because the X-set is a digital receiver and 
the 1 'Tsec, l/L’O bit samples are readily available from the carrier 
channel output. 
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Section 11. Data-Blt Synchronization 

After Costae lock has been achieved the receiver synchronises 
the user-time clock to the incoeiing data bits. The NOLL is tracking 
the C/A code which is synchronised to the data-bit transitions. How- 
ever the C/A code repeats every one millisecond whereas data-bit tran- 
sitions occur less frequently* every t%renty siilliseconds. Thus using 
only information from the NOLL, the receiver can accurately locate the 
data-bit transitions within a one-millisecond interval but is unable 
to determine which of twenty possible one -millisecond intervals in 
which the data-bit transition occurs. A histogram of the sum of the 
sign reversals of the inphase samples of the received signal in each 
of twenty one-millisecond time intervals is used to resolve this am- 
biguity. One second of inphase samples (or equivalently fifty data 
bits) is used to form the histogram. The sum of the sign reversals 
in the one-millisecond interval which coincides with the data-bit tran- 
sition must be larger than the sum in any other one-millisecond inter- 
val by at least ten sign reversals. When this is true, data-bit syn- 
chronization is declared. When data-bit synchronization has been de- 
clare) the user-time clock is appropriately advanced or delayed, to 
coincide with the data-bit transition. Then the integrate-and-dump 
interval is set to four milliseconds. The receiver then enters the 
data-frame synchronization mode. 

If at any time during the data-bit synchronization process Costas 
lock is lost, the procedure starts over. This prevents cycle skipping 
from being interpreted as data-bit transitions. 

A flow diagram of the data-bit synchronization process is shorn 
in Figure 24. 
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Secti on 1 2^ Datd">Frame Synchronisation 


After data*bit synchronization the receiver must be synchronized 

to the data frame being transmitted. A flow diagram of the data-frame 

synchronization process is snown in Figure 25. The process controller 

demodulates the data bits and compares them with the synchronization 

1 

bits (preamble) of the telemetry word and their logical inverse. (The 
telemc»try word is uriiformily distributed in the data frame and occurs 
every six seconds) . When the preamble or its logical inverse is detected 
the remaining bits in the telemetry word are demodulated and a parity 
check is performed. If a parity error is detected, the process starts 
over. If not, the handover word is demodulated and a parity check is 
made on it. As before, if a parity error is detected the process re- 
starts. If not, the current subframe is identified, the pseudo-range 
and user-time clock arc set, and tentative data-frame synchronization 
is declared. The remaining data until the next handover word is de- 
modulated, and a parity check is performed. If the Z-count (system 
time) of tl*e next handover word differs from the previous 2-count by 
more than one, data-frame synchronization is cancelled and the procedure 
is restarted. If the difference is exactly one, data-bit and data-frame 
synchronization are performed for the other channels that are Costas 
locked. Then the pseudo-ranges of the channels are compared. If any 
of them differ by more than 21 ms, data-frame synchronization is can- 
celled and the process is restarted. However, if the differences are 
less than 21 n;s , data-frame synchronization is declared. 

The conparison of pseudo-ranges is a reasonability test. Twenty- 
one milliseconds of the C/A code corresponds to approximately 4000 
miles. The maximum pseudo-range difference will exist when one of the 
satellites is directly overhead and one near the horizon. There is 
approximately 4000 miles pseudo-range difference between a satellite 
directly overhead and one 15® above the horizon. Normally satellites 
With elevation angles lower than 15® above the horizon will not be 
ust'd for navigation. Thus, if the pseudo-range difference is greater 
than 4000 miles the user-tim.e clock has probably been incorrectly set 
and tlie data-frame synchronization piocess (which sets the user-time 
r lock) IS repealed. 
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INTEROFFICE MEMORANDUM 
GPS-B1-CRC-3 

April 8. 1981 

LinCom Corporation From c. m. chie 

To; Dr. Jim Pawlowski 

Subject: Considerations for Testing GPS Receivers Using a Single GPS Signal Simulator 
cc; John MacLeod , Jack i^hnson 

Sunmary 

The attached package of “jtes summarizes a preliminary 
investigation of the problems and options associated with testing the 
GPS receiver. In particular, the emphasis Is on testing the GPS 
receiver when only one simulated emitter signal (out of a total of a 
complete set of four) is available. The first part of the package is a 
general description of the ideal test configuration and the basic system 
and subsystem tests required to determine receiver performance. 

The implications of testing with a nonideal configuration, for 
example, with 1 satellite simulator, are then considered. Special 
attention is devoted to the 2-channel sequential set (present JSC 
baseline). Supplementary tests that can be perfon;ied independent of the 
simulator signals are also highlighted. 

Additional Coriiients 

The material was discussed over the phone with Or. Jin Pawlowski on 
March Zb, 1981. At that titne several additional comments were made. 

>)c following is a sumnary. (Items 1 and 2 were discussed v/hile going 
i'lrough v/ays to fool a one channel sequential set.) 




(1) In order to reduce the Interface requirement of the R/PA with the 
outside world, it was suggested that the simulator be programmed 
to switch emitter signals independently at approximately the same 




I 

I 

t 


time as the R/PA. Figure 1 demonstrates the potential 
difficulties. 

(2) Another scheme for synchronizing the R/PA and the simulator is 
through an external timing source, i.e. , 



This requires the R/PA to accept external sequencing controls. 

(3) The R/PA design may want to incorporate some st.t of "data" 
testitig capabilities (i.e., provide test access ;>oints to the 
right of the dashed line in the R/PA functional block diagram in 
the package.) T.ie advantage is that the digital/conputational/ 
control portion of the R/PA can be checked out more rapidly and 
independent of the RF section. 

(4) Consider R/PA procedures for entering into, staying in and 
recovering fror! a blackout. 

(j) Consider R/i’A procedures when less than 4 emitters are visible. 

(o) t’ryvide capability for the R/PA to track a single satellite 
continuously for testing ivirposo. 

'. /) Consider P /PA test access points and externally controllable 


functions based on selected testing philosophy in the next round 
of spec revision. 
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« TYPES OF RECEIVERS 


• ] CHANNEL (SEQUENTIAL) 


• 2 CHANNEL (SEQUENTIAL) 


• 5 CHANNELS (PARALLEL) 


• SATELLITE SIMULATOR AVAILABILITY 


• 4 SIMULATORS 


• 1 SIMULATOR 


I SIMULATOR/RECEIVER INTERFACE REQUIREMENTS 


• DEPENDS ON ABOVE 


• nW.T ON RECEIVER DESIGN W.R.T. TEST ACCESS 
POINTS AND RECEIVER CONTROL 


• "FIELD" TESTS 


• TESTS riRFORMED r.'ITH ACTUAL SATELLITE SIGNALS 
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FOR COMPLETE TfSTING OF R/PA 


SYSTEM LFVFL TESTS 

• SATELLITE SELECTION 

0 TIME-TO-FIRST-FIX 

0 R/PA ACCURACY 

0 PSEUDO-RANGE 
0 RANGE-RATE 

0 NAVIGATION ACCURACY 

SUBSYSTEM LEVEL TESTS 

0 RECEIVER HARDWARE 
0 CODE LOOP 

0 ACQUISITION 
0 TRACKING 
0 CARRIER LOOP 

0 ACQUISITION 
0 TRACKING 
0 BIT SYNCHRONIZER 
0 FRAME SYNCHRONIZER 
0 CODERS 

0 USER TIME CLOCK 

0 RECEIVER CONiROL 
0 ACQUISITION 
0 CHANNEL SEQUENCING 

0 FREQUENCY PREPOSITIONING 
0 CODE PREPOSITIONING 
0 CONSTELLATION REVISION 

0 NAVIGATION SOFTWARE 
0 POSITION FIXES 
0 TIME FIXES 



WITH ONE SIMULATOR AVAILABLE 


• ONE CHANNEL SEQUENTIAL 

• FULL-SCALE TESTING POSSIBLE PROVIDED 
INTERFACE REQUIREMENTS ARE MET 

• TWO-CHANNEL SEQUENTIAL 

• PARTIAL TESTING ONLY 

t FULL-SCALE TESTING POSSIBLE WITH 2 SIMULATORS 

• 5 CHANNEL PARALLEL 


• PARTIAL TESTING ONLY 





SEQUENTIAL ONE-CHANNEL SET INTERFACE REQUIREMENTS 




• T IS THE TRANSITION TIME REQUIRED TO PREPOSITION THE 
FREQUENCY AND CODE FOR A SECOND EMITTER 

• S IS THE SWITCHING TIME REQUIRED TO CHANGE FROM ONE 
EMITTER SIGNAL TO ANOTHER 

• D IS THE TIME DELAY REQUIRED TO RELAY THE SWITCHING OF R/PA 
TO SIMULATOR 


• IF D+S < T, THE R/PA IS COMPLETELY FOOLED BY ONE SIMULATOR 








OPTION #1 


• FORCES 2 CHANNEL R/PA TO WORK AS TWO SEQUENTIAL 1 CHANNEL 
RECEIVER OPERATING ON L^ AND L^ 

CH 1 (L^) I EMITTER #1 | #2 | #3 | #4 | 

CH 2 (L^) I EMITTER #1 | H2 | <>3 I #4 I #1 

• COMPLETE TESTING POSSIBLE 





OPTION 02 
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• PROVIDE CONTROL CAPABILITY TO R/PA SO IT CAN TRACK SAME 
EMITTER MORE THAN ONCE 


CH 1 

SAT #1 

SAT n 

SAT #1 

SAT #2 

CH 2 

#1 

n 

#1 

#2 


• IMPACTS: 


• POOR SYSTEM LEVEL TESTING 

• POOR GDOP 

t RECEIVER HARDWARE TEST OK 

• RECEIVER CONTROL OK EXCEPT FOR CONSTELLATION REVISION 
t NAVIGATION SOFTWARE 

0 POOR GDOP 

















t ADDITIONAL TESTS NEEDED (WITHOUT USING SIMULATOR) 

• COMPUTER PROCESSING PORTION OF R/PA CAN BE TESTED 
BY COMPUTER SIMULATION (SEE FIGURE) 

• OVERALL SYSTEM LEVEL TEST MUST BE VERIFIED WITH 


ACTUAL GPS SIGNALS 
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Sumnary • 

In LInCom's last report [1] Me described a method of obtaining bit and 
subframe synchronliatlon. In this report we present the analysis and 
the performance of the bit synchronization technique presented to 
portray relationships between the SNR and the probability of correct 
detection of the bit synchronization circuit along with the average wait 
time for each decision. 

The results are presented for the relevant CNR ra.ige of 29 ^ C/Nq ^ 
36 dB*Hz which when converted to the bit signal to noise ratio SNR 
becomes 12 _< SNR _< 19 (fi. From the curves presented In this report one 
can say that for a probability of correct detector of 0.99 and an SNR of 

10.2 dB It would be necessary to have the shift register length of at 
least 5 and the average waiting time to make a decision would be about 

1.2 seconds. 

1. GPS Signal Structure and Bit Synch Circuit 

The Incoming GPS signal at any receWer has the C/A (clear 
acquisition code), the P (precision timing) code and the data message. 

One Important thing to notice here that the C/A code epochs (1000 epochs 
per second), P code epochs (1 epoch In every l.S seconds) are coherent 
with data bit boundary, (50 boundaries per second). Thus each data bit 
contains exactly 20 C/A epochs, I.e. within every count of 20 of the C/A 
epochs, one and only one C/A epoch Is coincident with the data bit 
boundary. In this report analysis will be carried on only for the bit 
synchronization circuit which decides which C/A epoch out of the 20 is 
coincident with the data bit boundary. 

The C/A code Is a pseudo random noise chip stream unique In pattern 
to each satellite and as was mentioned earlier repeats once every 
millisecond. It Is relatively easy for the receiver to match and lock 
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ACCUMULATOR I 
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onto the C/A code because the code search Is limited to time interval of 

f 

1 millisecond and the chip rate is only one tenth of that of the P 
code. Thus for the analysis to follow we assume that the local C/A code 
generator is in lock with the incoming code. Consequently, the C/A code 
epochs are at our disposal. 

Figure 1 shows the timing structure of the received GPS signal and 
Figure 2 shows the C/A bit synchronization resolution circuit. The 
working of the circuit can be described briefly as follows. 

The bit sync circuit operates by checking for data transitions 
around each epoch (20 (Kjssible epochs per bit). There are 20 
accumulators keeping track of the transitions scored by a particular 
epoch. After one of the accumulators fills up to 2® * 64 counts, the 
associated epoch is declared to be the data bit boundary and the data 
clock is set accordingly. The data clock is used to provide the timing 
for the 20 millisecond integrate and dump circuit for the data detector. 

In order to simplify hardware, the possible data transition is 
determined by looking at portions of 1/20 bit as shown. The transition 
detector shown in Figure 2 can be a simple comparator. Notice that 
since only 1/20 of a bit is integrated, the signal-to-noise ratio 
degrades by 13 <B. To compensate for this degradation, the accumulators 
are used. The signal-to-noise ratio improvement is proportional to the 
number of transitions observed before declaring "lock." 

2. Mathematical Modeling 

As said above, the transition detector has twenty different output 
ports, each output is recorded by an accumulator attached to it. We 
also know that for every bit time there is only one C/A epoch coinciding 
with the bit boundary. And if a transition in bits is present then that 
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port will put out a 1 (transition present) with a high probability ,(we 
will cjII this port the correct port) end at the same time all the 
remaining ports will have zeros with high probability. The procedure 
will repeat everytlme there Is a transition In bit putting a *1* In the 
accumulator attached to the correct port. This continues until the 
accumulator overflows announcing the boundary of the bit at the C/A epoh 
associated with the correct port. Assuming uncorrelatedness of the one 
bit decition and the transition detector decision frwn the current bit 
slice under observation to the previous bit slice it Is easy to write 
down the probability of a correct decision. Assuming '£* is the length 
(number) of I's at which the accumulator will overflow (ldg2 * * number 
of stages of the accumulator) 


1 


• £-1 

ProbCCorrect Decision] ■ I Jf I PS(l-Po)'‘’''(!iJ]' I 

k=Z m*0 


P, = Probability of detecting a transiton at the 
correct port 


pQ ® Probability of detection a transition at any 
one of the incorrect ports 


Assuming that transition and nontransition of bits are equally likely, 


we get 


Pj « 0.5 


where 


Pq - 2P(e)[l-P(e)3 

P(e) * Probability of deciding +1 given Tl was 
transmitted 


; _L 


Ejj/Nq Is the bit signal to noise ratio. 
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Next, we would like to find the average waiting time given that a 


correct decision is made. This can be written down as 

Average Waiting Time * I kpj (1-Pi J 

k=l 

= ~ (see Appendix) 

Pi 

this is a general result, but in our case we have pj = J.5; hence, the 
average waiting time for the shift register length of log£ t is 2t» 

Figure 3 shows the variation of P(C) with respect to threshold 't', 
SNR being the parameter. As can be seen, for low SNR (Fjj/Nq) such as 8 
dB even the threshold of 64 (i.e. shift register length of 6) gives 
probability of correct detection to be no more than 0.975. But for high 
SNR's say 15 dB, even a threshold of 16 gives a probability of correct 
detection to oe in excess of 0.9999. Figure 4 depicts P(C) vs SNR with 
the threshold £ to be the parameter. 

Figure 5 plots ' £* vs the SNR (dB) with P(e) to be the parameter. 
For example with the SNR of 11 dB and a desirable probability of correct 
detection of 0.9999 we can see that the threshold requirement is about 
42. 

It should be pointed out that throughout this analysis the 
uncorrelatedness of the decision making process for the adjacent bit 
slices was assumed resulting in uncorrelatedness of the consecutive 
input of the comparator. Actually there does exist some nonzero 
correlation between the two consecutive outputs of the correlator which 
makes our results slightly pessimistic. 
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3. Cone Iws ions 

For a probability of correct detection of the bit boundary to be 
0.99 and a reasonable F5/N0 of about 10 dB it will require a shift 
register length of the accumulators to be 5. This arrangement will 
require on the average 64 data bits or about 1.2 seconds of time to 
detect the data bit boundary. 

At j irw SNR such as 8 dB to obtain the probability of correct 
detection to be 0.99 will require 7 shift registers in the accumulator 
and which in term imposes an average waiting time of about ?56 data bits 
cr about 5 seconds. 



APPENDIX 


Let 


S 


t 


1 




k-1 

£-1 


J 


It can be shown that 


vk-£(k-l 


Pi ,£ r 


s , ■ ,i p ,‘( i - p , rr.:jj = Piiir^r ^( i-pi 


k*£ 


Pi k=l 


let 1-pi = qj 


s, = Pi(^j' I 

» ^ Ml l. = £ ' ‘ 


Pi V k 
— 1 qi k ’ 

‘’i k=i " 

? d . 1 I 
Pi dq ‘T-qj ' 

Pi ,£♦! 


£L Ip JL y nh 

qi <^1 dq, J , ^ll 

Pi 




p^ 4 


p'l,£4l £+1 ^ k, k ,, fPlit^lj V V ^ I 

= Pita-J q ^ I qiii+lJ ' Pl^^J I ^ ‘I ^ ^ ^£* 1 '' 

^ ^1 k-£+2 ^ ‘ ^ ‘’l k»£ " k = £ + l 


k + lf k 


’ * q,Sj„ 


giving us 



Since qj = 1-pj we have 




Pi ,t+i 


= 1 


k=« 


Now we will proceed ahead to find • 


Average 
Waiti ng 
Time 


k* 1 


we know that 

i Pi'(i-Pi)'^‘‘ij;ij ■ 1 

k=£ 





0 


i ‘p,‘''(i-p,)(i:Jj - i pf(k-i)(i-p)'‘‘'' 

k-£ ^ k'£ 


= 0 


«(i * 1^) Pi‘0-Pi)'''‘lil!) ' 11^5 Ji ‘‘t-l) 





Therefore 


[1] c 
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